A unified dataset for the city-scale traffic assignment model in 20 U.S. cities

City-scale traffic data, such as traffic flow, speed, and density on every road segment, are the foundation of modern urban research. However, accessing such data on a city scale is challenging due to the limited number of sensors and privacy concerns. Consequently, most of the existing traffic datasets are typically limited to small, specific urban areas with incomplete data types, hindering the research in urban studies, such as transportation, environment, and energy fields. It still lacks a city-scale traffic dataset with comprehensive data types and satisfactory quality that can be publicly available across cities. To address this issue, we propose a unified approach for producing city-scale traffic data using the classic traffic assignment model in transportation studies. Specifically, the inputs of our approach are sourced from open public databases, including road networks, traffic demand, and travel time. Then the approach outputs comprehensive and validated citywide traffic data on the entire road network. In this study, we apply the proposed approach to 20 cities in the United States, achieving an average correlation coefficient of 0.79 in average travel time and an average relative error of 5.16% and 10.47% in average travel speed when compared with the real-world data.

which is a mature model that has been studied extensively in the transportation field, is adopted to estimate the city-scale traffic states.The input of the traffic assignment model only includes the Origin-Destination (OD) demand information and network structure, both of which are public and openly available.Then, the model outputs the city-scale traffic dataset.Traffic assignment models utilize OD data to predict traffic flow and route choices for individual travelers, relying on either predefined or data-driven behavioral models.By modeling the interactions between travelers' behaviors and traffic congestion, the traffic assignment model searches for the equilibrium condition that mimics real-world traffic conditions.Traffic assignment models can often serve as the primary tool for local governments to assess the potential impact of changes in land use or road network expansions on both local and global traffic conditions.These models are indispensable because they inherently focus on optimizing travel decisions for local residents, aligning with their individual preferences.This capability enables the model to predict changes in agent-level behavior in situations that may not be fully reflected in the available data.Moreover, traffic assignment models demonstrate robust predictive capabilities for estimating future traffic conditions.For example, Metropolitan Planning Organizations (MPOs) in urban areas of the United States would utilize travel survey data, such as the National Household Travel Survey (NHTS: https://nhts.ornl.gov/), to produce traffic data for each local urban area that represent residents' travel patterns 19 .However, these traffic assignment models and data are usually maintained by public agencies and generally not available to most researchers or the public due to difficulties in information sharing or privacy concerns 20,21 .Furthermore, the data used in traffic assignment models are under the ownership of various institutions and lack standardization in terms of their structures, granularity, and output formats.As a result, the data are restricted to a few researchers and it is challenging to access the necessary data for traffic assignment models across cities from official sources.Given the above, there is still a notable absence of city-scale traffic datasets that include multiple major cities within one geographic and cultural region, adhere to consistent standards, collect and validate information on a uniform scale, provide comprehensive data types, and meet high-quality standards for public availability.
Although there are a few publicly available datasets 22,23 concerning urban areas (see Table 1), the reliability and completeness of these datasets limit their applications across broader urban studies, especially in fields like energy, environment, and public health 24,25 .The limitations come from the following aspects: First, the existing traffic datasets typically cover some important traffic segments for a single city rather than a city-scale traffic dataset for multiple cities.Second, these current datasets often lack the necessary input, including road network data and corresponding OD data, directly usable for traffic assignment models.Third, these datasets often suffer from incomplete data types and lack of timely updating, resulting in limited convenience when utilizing them.In other words, these datasets are often collected by different researchers or volunteers several years ago, leading to a lack of uniformity in the data types and formats, as well as infrequent updates and maintenance.Fourth, these datasets frequently lack comprehensive validation across multiple variables or fail to offer adequate tools for predicting traffic features from behavioral data.For example, a dataset that includes OD numbers may result in unrealistic traffic flow predictions when attempting to utilize a traffic assignment model.In light of these mentioned facts, currently, there is no unified and well-validated traffic dataset available for multiple cities that covers the entire urban road network at a citywide scale, which hinders the feasibility of conducting comprehensive urban studies across cities to unearth novel discoveries.
To facilitate convenient access to citywide traffic assignment models and data for researchers from different domains besides transportation fields, this study provides a unified traffic dataset for traffic assignment models in 20 representative U.S. cities, with populations ranging from 0.3 million to over 8.8 million.Specifically, we first obtain the input of the model by fusing multiple open public data sources, including OpenStreetMap, The Longitudinal Employer-Household Dynamics Origin-Destination Employment Statistics (LODES), Waze, and TomTom.Then, we employ a grid-search method to fine-tune the parameters and generate the final traffic dataset for each city.The real world's average travel time and traffic speed serve as validation criteria to ensure a reliable and effective traffic dataset for multiple cities.The validation results demonstrate that our approach can successfully produce the dataset with an average correlation coefficient of 0.79 for average travel time and an average error of 5.16% and 10.47% for average travel speed between real-world data and our data.Finally, we upload the validated traffic dataset and the code used in this study to a public repository.To sum up, we utilize the static traffic assignment model, leveraging annually aggregated statistical data and open public data sources, to offer a city-scale traffic dataset for macroscopic urban research.It is worth noting that the approach provided in this study can also be applied to other cities.A comprehensive workflow of processing multi-source open public datasets to acquire this dataset is provided in Fig. 1.

Methods
Creating a unified traffic dataset in multiple cities involves four main procedures: (1) the identification of representative cities; (2) the acquisition of corresponding input data from multi-source open public datasets; (3) the fusion of the obtained data; and (4) the implementation of traffic assignment, along with parameters calibration.The main procedures are illustrated accordingly below.

Identification of representative cities.
In this study, we select a total of 20 representative cities in the United States and generate corresponding traffic datasets using the proposed approach.To ensure diversity and exemplarity among the selected cities, we primarily consider factors such as geographic location, urban scale, topography, and traffic conditions during the commute.Our selection includes a range of cities, including megacities like New York City, as well as several large cities such as Chicago and Philadelphia.We also included smaller but equally representative cities such as Honolulu.The topography of these cities also varies widely.For example, New York and San Francisco are separated by several rivers and rely on critical bridges and tunnels for commuting, while Las Vegas and Phoenix have relatively flat and continuous terrain, with surface transportation playing a predominant role.
Basic information of the 20 representative cities in the United States is given in Table 2.The population and land area data in the year 2020 are sourced from the U.S. Census Bureau (https://www.census.gov/)while the congestion ranking information in the year 2022 is from TomTom (https://www.tomtom.com/traffic-index/ranking/).Their geospatial distribution is shown in Fig. 2.    crowd-sourced road network geographic information, including network topology, road attributes, and connectivity information.By leveraging OSM data, researchers gain convenient access to a comprehensive and up-to-date depiction of the network structure, which facilitates the research in urban studies [26][27][28][29] .The road attributes are also sourced from OSM.After the implementation of cleaning and integration procedures, these processed data can serve as the input for the traffic assignment.A summary of the road network data is given in Table 3. Specifically, we employ a Python package named osmnx 30 (https://github.com/gboeing/osmnx) to download the OSM data.We then use another Python package called osm2gmns 31 (https://github.com/jiawlu/OSM2GMNS) to extract the nodes and links on the road network from the OSM data and save them into separate CSV files in GMNS format 32,33 .We use five main link types including 'motorway' , 'trunk' , 'primary' , 'secondary' , and 'tertiary' to implement the traffic assignment.For each link type, we initiate the corresponding road attributes, including parameters such as road capacity, speed limits, the number of lanes, and so on.For the nodes, each node represents the intersection between two links and contains a unique identifier along with latitude and longitude information.By establishing the connectivity between nodes and links through their corresponding relationships, the network topology and road attributes can be constructed.We use the graphing functions of osmnx to visualize the constructed road networks of 20 representative U.S. cities in Fig. 3.
Travel demand.We then estimate the travel demand, another essential input data for traffic assignment, using the data from the LODES dataset (https://lehd.ces.census.gov/data/lodes/)provided by the U.S. Census Bureau.The LODES dataset includes commuting data for the workforce in all states across the United States over multiple years, which have been widely used in existing studies 34 .LODES data collection involves employers reporting employee details to state workforce agencies, including work and home locations.The U.S. Census Bureau collaborates with state agencies to process and anonymize this data.It's then used to create Origin-Destination (OD) pairs.This dataset, at the finest granularity of block level, documents the block code for both workplace census and residence census, along with the corresponding total number of jobs.Essentially, the LODES dataset provides an excellent representation of the trip distributions of the U.S. working population that can be used to construct the OD matrix.In this study, we mainly focus on producing the traffic dataset for the year 2019 and the commuting OD data in that year are collected.Moreover, the data collection process is performed at the block level, resulting in the OD data between blocks.
Travel time and speed.We collect data from two open-source dataset platforms, namely TomTom (https:// www.tomtom.com/traffic-index/ranking/)and Waze (https://www.waze.com/live-map/),as two indicators of travel time and average speed respectively for our dataset validation.The detailed procedures of data collection can be found in the subsequent sec:Technical ValidationTechnical Validation section.

Data fusion.
In this section, we integrate the road network data and OD data to unify the data format.Since the origins and destinations in the OD matrix are not associated with network nodes, it is infeasible to directly take these data as input for the traffic assignment.Therefore, we need to establish a connection between network nodes and blocks.After establishing the connection, we can employ the traffic assignment model to identify appropriate travel paths and allocate traffic flow to the respective links.To be specific, we begin by aggregating the OD data from its minimum granularity at the block level to a higher level, namely, the tract level.According to the United State Bureau [35][36][37] , blocks are statistical units with small areas, generally defined to contain between 600 and 3,000 people, whereas tracts composed of multiple blocks are relatively larger and typically have a population size ranging from 1,200 to 8,000 people.In order to achieve a balance between computational complexity and accuracy, we consider tracts as an ideal basic unit for the traffic assignment, which is similar to the existing studies 38,39 .This implies that we use the tract as a Traffic Analysis Zone (TAZ) in the traffic assignment model.
Then, the geographical location of each TAZ is determined as the average coordinates of all the blocks within a tract.These TAZs (also called centroids) are generated and stored in the existing node file labeled with a unique identifier.Finally, we generate connectors to bridge the TAZs and network nodes.These connectors can be regarded as a special type of links that are generated from each TAZ center to their neighbor links.Moreover, these connectors are incorporated into the existing links labeled with a unique identifier.As a result, the commuting trips could start from the origin TAZ, traverse a connector to access the nearby road network, choose a suitable path, and then use another connector to reach the destination TAZ.
Traffic assignment.In this section, we use the traffic assignment model to produce the dataset based on the User Equilibrium (UE) 40 .To be specific, we formulate the UE using an optimization model and calibrate four categories of parameters used in the model.Using the network structure and OD demand as input, the model would output the traffic flow, speed, and density on each link.Moreover, we mainly focus on the static traffic assignment and do not consider the influence of temporal variations on traffic conditions.User equilibrium.All travelers naturally make decisions to minimize their own travel costs (either travel time or equivalent monetary value).Wardrop's First Principle 41 posits that when every traveler seeks to minimize their individual travel costs, traffic flow eventually stabilizes.In this equilibrium state, the travel costs on all utilized paths become equal and minimized.Meanwhile, the travel costs on unused paths for any given OD pair are greater than or equal to those on the used paths.In other words, a steady-state traffic condition is reached only when no traveler can improve his or her travel time by unilaterally changing routes.The satisfaction of Wardrop's first principle is commonly referred to as User Equilibrium (UE).
The physical transport network including road segments and intersections in an urban area can be represented as a graph structure G(N, A) containing a link set A and a node set N. For each link α ∈ A, it has the link flow x a and the link travel cost t a respectively.For each node r, s ∈ N, it is defined as the TAZ that generates or attracts traffic demand.Therefore, the mathematical formulation of the traffic assignment model under the UE condition 42 can be expressed as follows: where t a (x a ) denotes the link performance function that indicates the travel cost on link a when the traffic flow is x a .f k rs represents the traffic flow on path k connecting origin r and destination s. q rs indicates the number of trips from origin r to destination s. δ ka rs is a binary variable indicates whether link a is part of path k between origin r and destination s.Equation ( 2) imposes the flow conservation constraints.Equation (3) expresses the relationship between link flow and path flow.Please refer to the book Urban Transportation Networks 40 for details.
Once the traffic flow on each link is determined, the total travel time, denoted as c k rs , for a specific path k can be calculated by summing the travel time of each link along this path, which can be formulated as follows: k rs a a ka rs Although the above optimization model has been proven to be a strict convex problem with a unique solution for traffic flow on links 40 , the computational cost of finding the optimal solution would significantly increase when dealing with large-scale city road networks.To alleviate the computational burden, a bi-conjugate Frank-Wolfe algorithm 43,44 is employed to find the optimal solution.In order to enable convenient usage of the provided dataset by users from various disciplines and allow them to easily modify the core parameter settings of the traffic assignment process according to their research needs, we employ two traffic modeling platforms to generate the final dataset.Subsequent users can either directly view the dataset in a no-code format or quickly adjust parameters through a low-code approach to conduct scenario testing under different scenarios.Specifically, a commercial software (named TransCAD) and an open-source Python package for transportation modeling (named AequilibraE) are utilized simultaneously in this study.For both platforms, the maximum assignment iteration time and the convergence criteria are set to 500 and 0.001, respectively.The results of the traffic assignment model in 20 U.S. cities are shown in Fig. 4.

Parameters calibration.
The traffic conditions on the network are influenced by many factors related to traffic supply and demand.Consequently, the traffic assignment model would be impacted and output different results.Since the disturbances in the transport system are nonlinear and challenging to quantify, it is difficult to establish a deterministic mapping relationship between various influencing factors and the results of the traffic assignment model.Therefore, we adopt a grid-search approach to calibrate four common categories of factors that are closely related to the traffic assignment model.We determine the final model by continuously fine-tuning various parameters associated with the traffic assignment model until the transport system reaches the UE condition.In this study, we introduce four categories of factors including road attributes, travel demand, impedance function, and turn penalty, as outlined below.

• Road attributes
We categorize the entire road network into three major types, namely expressways, arterial highways, and local roads.Capacity and free flow speed of each road type are two parameters identified to be calibrated.
Based on the experimental results, the appropriate range of road capacity for expressways is between 1800 veh/h/lane and 2200 veh/h/lane, while the range for free flow speed is from 65 km/h to 90 km/h.In the case of highways, the corresponding capacity value falls within the range of 1500 veh/h/lane to 2000 veh/h/ lane, and the free flow speed value ranges from 40 km/h to 65 km/h.As for local roads, their capacity varies from 600 veh/h/lane to 1500 veh/h/lane, while the suitable speed ranges between 25 km/h and 45 km/h.The detailed information for each type of road can be found in Table 4.

• Travel demand
The OD travel demand is another significant factor influencing the outcome of the traffic assignment.In this study, we aim to simulate the traffic conditions during the peak hours.As mentioned above, the OD demand matrix is derived from the total number of jobs in the United States in 2019, generated from LODES datasets.Although it is reasonable to assume that commuting travel accounts for the majority during peak hours, such demand cannot reflect the actual traffic conditions.Therefore, it is necessary to adjust the initial OD demand, considering variations in transport modes, travel departure time, and carpooling availability during commuting to work.To address this issue, we introduce an OD multiplier to estimate the actual traffic demand during the commuting time.We find that stable results can be obtained when the parameter ranges from 0.55 to 0.65.We show the travel demand and the percentage of internal travel within each TAZ in Fig. 5.

• Link performance function
The link performance function, also known as the impedance function or volume delay function, refers to the relationship between travel time and traffic flow on a road.Typically, travel time increases non-linearly with the increase in traffic flow, which also significantly affects the traffic assignment.One of the most commonly adopted functions in the literature is called the Bureau of Public Roads (BPR) function 45 , which is expressed as follows: are parameters needed to be fine-tuned.We find that the results are satisfactory when parameter α ranges from 0.15 to 0.6 while parameter β changes from 1.2 to 3. The specific values of parameters for each city are provided in Table 5.

• Turn penalty
The turning delay at intersections is also a significant factor that should not be dismissed.When vehicles pass through road intersections, their speed typically decreases, either due to signal control or the necessity to make turns.However, this behaviour cannot be adequately represented in solving traffic assignment problems.To ensure that the results of the traffic assignment model are in accordance with real-world scenarios, we uniformly set corresponding parameters for all junctions to simulate the turning delay effects.In other words, the turn penalty parameters are an average value for the turning delay at all intersections in the road network and these intersection types include signal-controlled intersections, roundabouts, yield or stop intersections, and others.Specifically, the time delay for right turns varies between 0 and 0.25 minutes, while the penalty for making a left turn ranges from 0 to 0.35 minutes.The delay for through traffic is between 0 and 0.15 minutes.U-turn is prohibited in the traffic assignment simulation.The specific parameter setting is demonstrated in Table 5.

Data records
We share the traffic dataset on a public repository (Figshare 46 ).In this dataset, each folder, named after the city, contains the input and output of the traffic assignment model specific to that city.We elaborate on the details as follows: Input data.This folder contains all the input data required for the traffic assignment model, namely the OD demand data and network data.The network data contains both node and link files in a CSV format.The data in this file folder specifically includes the following contents: • the initial network data obtained from OSM • the visualization of the OSM data • processed node/link/od data The detailed meanings of the fields contained in different input data are given in Table 6.
TransCAD results.This folder contains all the input data required for the traffic assignment model in TransCAD, as well as the corresponding output data.The data in this file folder specifically includes the following contents: • cityname.dbd:geographical network database of the city supported by TransCAD • cityname_link.shp/cityname_node.shp:network data supported by the GIS software, which can be imported into TransCAD manually • od.mtx: OD matrix supported by TransCAD • LinkFlows.bin/LinkFlows.csv:results of the traffic assignment model by TransCAD • ShortestPath.mtx/ue_travel_time.csv: the travel time (in minutes) between OD pairs by TransCAD The detailed meanings of the fields contained in output data generated from TransCAD are given in Table 7.
Aequilibrae results.This folder contains all the input data required for the traffic assignment model in AequilibraE, as well as the corresponding output data.The data in this file folder specifically includes the following contents: • cityname.shp:shapefile network data of the city support by QGIS or other GIS software • od_demand.aem:OD matrix supported by AequilibraE • network.csv: the network file used for traffic assignment in AequilibraE • assignment_result.csv:results of the traffic assignment model by AequilibraE The detailed meanings of the fields contained in output data generated from AequilibraE are given in Table 8.

Technical Validation
To ensure the consistency between the traffic assignment model's output and real-world traffic conditions, we conduct validation using two different public open sources of traffic data.Specifically, the travel time between different OD pairs and the overall average travel speed are employed as two validation indicators to ensure the reliability and accuracy of the provided dataset.The validation results are shown in Tables 9, 10 and we can see that the provided dataset for each city is accurate and valid.

Travel time.
In examining the travel time metric, we obtain the travel time between different OD pairs both from traffic assignment models and map service providers.As for the model side, the travel time under both UE and free flow conditions are calculated respectively using traffic assignment models.First, under UE conditions, the travel time between different OD pairs could be generated by summing the link travel time determined by the corresponding assigned traffic flow along the shortest path as shown in Eq. ( 5).Then, under free flow conditions, the travel time between OD pairs is the travel time associated with the shortest path, disregarding congestion on road segments.Furthermore, the average value of Travel Time (in minutes) under UE conditions (UETT) as well as free flow conditions (FFTT) for all OD pairs can be expressed as follows: where c ue rs and c ff rs denote the travel time between origin r and destination s under the UE and free flow conditions respectively.Additionally, the difference as well as the ratio between these two types of travel time give the average travel delay (in minutes) and delay factor for each city.
In terms of the real-world data for validation, since nowadays many map service providers have the capability to offer travel time estimates between two location points at different departure times based on users' historical navigation records, in this study, we choose Waze as the data source to obtain the actual travel time between each OD pair by using its WazeRouteCalculator API (https://github.com/kovacsbalu/WazeRouteCalculator)with Python code.
The results of travel time are shown in Table 9.It can be seen that Honolulu experiences the least travel time under free flow conditions, at about 8.70 minutes, while Minneapolis has the shortest average travel time during commuting hours, at about 10.25 minutes.Minneapolis also has the lowest delay travel time among all cities, merely 0.47 minutes, indicating that the commuting travel time in this city is almost the same as the travel

Data
Field Description

Output link_id
The unique identifier generated automatically by AequilibraE.matrix_ab The output assigned traffic flow (veh/h) for the AB direction.matrix_ba The output assigned traffic flow (veh/h) for the BA direction.matrix_tot The total output assigned traffic flow (veh/h) for both directions.

Congested_Time_AB
The congested link travel time for the AB direction.

Congested_Time_BA
The congested link travel time for the BA direction.
Congested_Time_Max The maximum congested link travel time for both directions.

Delay_factor_AB
The ratio of congested travel time to free flow travel time for the AB direction.

Delay_factor_BA
The ratio of congested travel time to free flow travel time for the BA direction.

Delay_factor_Max
The maximum ratio of congested travel time to free flow travel time for both directions.

VOC_AB
The Volume to Capacity (V/C) ratio for the AB direction.

VOC_BA
The Volume to Capacity (V/C) ratio for the BA direction.

VOC_max
The maximum V/C ratio for both directions.

PCE_AB
The output assigned traffic flow (PCE/h) for the AB direction.

PCE_BA
The output assigned traffic flow (PCE/h) for the BA direction.

PCE_tot
The total output assigned traffic flow (PCE/h) for both directions.time under free flow conditions.In contrast, New York City experiences significant delays, with a delay time of 24.47 minutes, revealing that the travel time during peak periods in New York is more than double that of the free flow condition.In terms of the delay factor, New York City has the highest value, reaching 2.24, followed by Chicago with a value of 1.65.Minneapolis and Pittsburgh have the lowest delay factor values, both at 1.05.
To evaluate the results, we use the Pearson Correlation Coefficient (PCC) 47 to measure the correlation between the actual travel time and the travel time produced by our model.The PCC r xy is defined as follows: where r xy denotes the Pearson's Correlation Coefficient.x i and y i are the individual sample points indexed with i. n represents the sample size.
Since the turning penalties are not incorporated in the traffic assignment algorithm of AequilibraE, the parameter settings in TransCAD and AequilibraE are not identical.Consequently, results of the two platforms are not entirely consistent.Considering the more comprehensive parameter settings in TransCAD, we adopt the results of TransCAD as the primary benchmark.We perform PCC analysis between Waze and TransCAD, as well as between TransCAD and AequilibraE, with the evaluation results presented in Table 9.
From the correlation analysis, we can find that all correlation coefficients R 2 are greater than 0.7, which confirms the accuracy and reliability of the results to some extent.We also visualize the correlation coefficient for each city in Fig. 6.It can be seen that the simulated travel time is consistent with the travel time in the real world.
Average speed.The overall average speed of the entire road network is another important indicator for validation.In this study, we use the speed data collected from TomTom Traffic Index as the actual speed to validate our model.We first calculate the average link-based speed of our model through dividing Vehicle Hours Travelled (VHT) by Vehicle Kilometers Travelled (VKMT).Then, the average OD-based speed values are derived from the ratio of distance to travel time between each OD pair.The Mean Absolute Percent Errors (MAPE) and Mean Absolute Errors (MAE) for both the link-based speed and the OD-based speed are used to measure the reliability of our model: where y i is the actual observed value, y i is the predicted value, and n is the number of samples.The results are summarized in Table 10.We find that the average MAPE and MAE values for the link-based speed metric are 5.16% and 1.77 km/h, respectively.Moreover, the average MAPE and MAE values for the OD-based speed indicator are 10.47% and 3.82 km/h, respectively.This implies that our approach can produce satisfactory and reliable results.

No. City
Network traffic impact on model performance.To validate the effectiveness and robustness of our model across cities, we further investigate how traffic conditions of a city affect the model performance.The MAE and MAPE values for link-based average speed metrics obtained in Table 10 are used to evaluate the model performance.The traffic conditions are characterized by two different indicators.One is the ratio of the total OD travel demand to the number of links for the entire road network, which can characterize the average OD demand and represent the traffic conditions of a city.The other is the average speed (km/h) in rush hour obtained from TomTom (refer to Table 10).If the values of average traffic demand are large, it reveals a congested city network experiencing substantial traffic demand, exemplified by cities like New York and San Francisco.Conversely, a small value suggests a city road network with low traffic demand, as observed in cities like Atlanta and Dallas.We can draw similar conclusions with respect to the average traffic speed.
The results are shown in Fig. 7.The red dashed line represents the linear regression trendline that has been fitted to the data points.The R 2 values of Fig. 7a and Fig. 7b are 0.0049 and 0.0218, respectively.This implies that there is no evident relationship between the model performance and the varying traffic demand of the network.Similarly, the R 2 values of Fig. 7c and Fig. 7d are 0.0212 and 0.0177, respectively.This suggests that the model performance is not affected by the varying traffic speeds in different cities.In summary, the proposed model exhibits low sensitivity to variations in city traffic conditions and achieves satisfactory performance across cities.

Usage Notes
The acquisition of OD data is crucial in performing the traffic assignment and producing the citywide traffic dataset.In this study, we utilize the commuting OD data (LODES) provided by the U.S. Census Bureau to generate the OD matrix.For cities in other countries, OD data can be substituted with alternative open data sources, such as OD data provided by TomTom (https://developer.tomtom.com/od-analysis/documentation/product-information/introduction).
Moreover, we use the average traffic time and average travel speed between different OD pairs in the real world to validate the results of our approach, ensuring its reliability and accuracy.If additional data sources are available, such as traffic flow data obtained from traffic detectors, we can also use the corresponding data to further evaluate the effectiveness of the provided dataset.
It is worth noting that the provided dataset is mainly used for macroscopic urban research and policy development across interdisciplinary studies.In view of this, the given dataset provides full spatial coverage of the entire road network, unlike existing traffic datasets that focus on specific areas.Hence, the provided traffic dataset and existing traffic datasets complement each other, which can better facilitate research in urban studies.Specifically, the full spatial coverage of the provided dataset makes it valuable for comprehensive macroscopic urban research and policy development, making a notable contribution to the literature, such as public transport planning, road expansions, the determination of bus routes, the estimation of the transport-related environmental impact and so on.In contrast, existing traffic datasets (e.g., PeMS) may exhibit incomplete spatial coverage, making them less suitable for the aforementioned macroscopic urban studies.Actually, the datasets containing fine-grained temporal information are more suitable for investigating regional traffic dynamics by leveraging the spatiotemporal relationship between the traffic data, such as traffic prediction, spatiotemporal propagation of shockwaves, calibration of fundamental diagrams, traffic data imputation, and so on.
In this study, the provided dataset lacks fine-grained temporal information due to the limited availability of input data.To fully understand dynamic traffic patterns, it is essential to consider both spatial and temporal dimensions within the traffic data.Consequently, developing a dynamic traffic assignment model that effectively captures the spatiotemporal interdependencies of traffic data is important.Moreover, employing daily traffic data for more fine-grained validation would enhance further urban research.

Fig. 1
Fig. 1 The workflow of obtaining unified and validated traffic datasets from multi-source open public datasets.

Fig. 5
Fig. 5 Total travel demand and the percentage of internal travel demand for 20 U.S. cities.

Fig. 7
Fig. 7 The model performance in relation to different traffic conditions for 20 U.S. cities.(a) The MAPE values (%) regarding the average OD demand for different cities.(b) The MAE values (km/h) regarding to the average OD demand for different cities. (c) The MAPE values (%) regarding the average speed for different cities.(d) The MAE values (km/h) regarding the average speed for different cities.

3 .
AequilibraE_assignmnet.py: A Python code file for users who want to get access to the traffic assignment results by using the AqeuilibraE.

Table 1 .
Comparison between existing open public traffic datasets and our dataset.

Table 2 .
Basic information on 20 representative U.S. cities.

Table 3 .
A summary of the road networks data for 20 U.S. cities.

Table 4 .
Capacity (veh/h/lane) and free flow speed (km/h) of roads.

Table 5 .
Parameters for BPR function and turn penalty (minutes).

Table 8 .
Fields description for the AequilibraE output data folder.

Table 9 .
Correlation analysis for the average travel time (min) among Waze, TransCAD and AequilibraE.

Table 10 .
Comparison of the average speed (km/h) generated from our model and TomTom Traffic Index.